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The activity of the soft gamma ray repeater SGR 1806-20 is studied in correlation with the 
EXPLORER and NAUTILUS data, during the year 2004, for gravitational wave (GW) short signal 
search. Corresponding to the most significant triggers, the bright outburst on October S**^ and 
the giant flare (GF) on December 27*'', the associated GW signature is searched. Two methods 
are employed for processing the data. With the average-modulus algorithm, the presence of short 
pulses with energy Eg^ > 1.8 ■ 10*^ erg is excluded with 90% probability, under the hypothesis of 
isotropic emission. This value is comparable to the upper limits obtained by LIGO regarding similar 
sources. Using the cross-correlation method, we find a discrepancy from the null-hypothesis of the 
order of 1%. This statistical excess is not sufficient to claim a systematic association between the 
gravitational and the electromagnetic radiations, because the estimated GW upper limits are yet 
several orders of magnitude far away from the theoretically predicted levels, at least three for the 
most powerful SGR fiare. 
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PACS numbers: 04.80.Nn, 95.55.Ym, 95.85.Sz 



I. INTRODUCTION 

Discovered since 1979 [212], despite their unpredictable 
behavior, the soft gamma repeaters (SGRs) are consid- 
ered interesting targets in gravitational wave (GW) stud- 
ies. The fundamental ideas are the neutron star (NS) 
nature of these objects and the association of SGRs with 
magnetars [21 |4]. Detected as persistent X-ray source 
at ~ lO^^ergs s~^, they occasionally emit energetic soft 
gamma bursts, up to ~ lO'^^ergs s~^, or even much more 
energetic events. The large number of observed char- 
acteristics of SGRs, including the bursting activity dur- 
ing the three giant flares (GFs) detected to date [T| 0- 
[7], confirm the NS nature of SGR, and offer an effec- 
tive evidence of the presence of very high magnetic field 
{B ~ lO^^G). (For a review, see [8]). According to several 
models, the burst triggers are primary reference studying 
the gravitational radiation [51-fT5]. The most optimistic 
GW emission model [TU] foresees an extreme bursting 
activity, probably due to a sudden global reconfiguration 
of the internal magnetic field, and the magnetar shape 
deformation, causing the increase in the momentum of 
inertia with a release of energy through gravitational ra- 
diation. Also the proximity (^ 10 kpc) makes these ob- 
jects intriguing for GW searches, enhancing chances of 
detectability. 

To this aim, several measurements had previously been 
performed. Studying the GF occurring on December 27*'* 
2004, the AURIGA group 14 explored the frequency 
range 930-935 Hz, under the hypothesis of oscillating 
emission with a damping time of 100 ms. Expressing the 
result in terms of the dimensionless amplitude h, they 
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found an upper limit of 2.7 • lO^^'^, at the time of the 
hyperflare. In relation to the same event, the LIGO col- 
laboration [T5| examined the pulsating tail of the burst 
which revealed the presence of quasi-periodic oscillations 
(QPOs) in the X-ray light curve, as RXTE and RHESSI 
satellite had detected [T6l[l7j. LIGO found no excess and, 
racing from ^ lO''^ ergs, they set several upper limit lev- 
els on the GW emitted energy, depending on time and 
frequency radiation. More recently, the same LIGO col- 
laboration presented the results of short-duration GW 
events associated with SGR 1806-20 and SGR 1900+14 
storms occurred during the year starting from November 
2005 [T5H2U] . Including the GF on December 27*'*, they 
analyzed almost two hundred events finding no evidence 
of any association. Depending on signal waveform types, 
several upper limit are estimated, referring to more than 
five orders of the GW energy magnitude (from 10^^ erg to 
almost 10^"'^ erg) for hypothetical isotropic GW emission. 

In the present paper we consider the particular phe- 
nomenon of the SGR 1806-20 during the year 2004. The 
resonant GW detectors EXPLORER and NAUTILUS 
[H] data are analyzed corresponding to the long sequence 
that occurred on October 5*'' [35], and to the GF on 
December 27*''. As previously seen, during this event, 
several observations were performed using GW experi- 
mental apparatuses but, as far as the October outburst 
is concerned, there were no previous GW measurements 
although it represents a very significant event because of 
the exceptional energy release and duration, extending 
3 • 10^^ erg over more than one hundred pulses in a ten 
minutes span of time. 

The principal points of the present analysis are as fol- 
low. The presence of short pulses (up to a few tens of mil- 
liseconds) was searched for on both GW detectors, using 
the data filtered by linear filter matched to S-like signals. 
In Sec|ll| the sensitivity of the apparatuses is shown rel- 
atively to the year 2004. As explained in Sec jlll A] the 
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physical regions were chosen on two long sequences of 
twenty-six plus forty-six seconds, corresponding to the 
highest EM peaks of the October 5*'* outburst, and on 
one second centered around the trigger time of the GF 
on December 27*''. The long period, in the case of the 
outburst, and knowledge of the arrival time on Earth 
of a light-speed signal at the time of the GF [71 [Tl] 
are circumstances particularly opportune to avoid un- 
certainty in the analysis focusing on the reference time 
of the gravitational radiation. The expected distribution 
under the null-hypothesis is built using real data, with 
random time shifts of one GW detector with respect the 
other one (Sec IIIB). Essentially, two algorithms are em- 
ployed. With the first, discussed in Sec|IV[ we evaluate 
the excess amplitude on the on-source regions of the GW 
detector data. With this measurement (calibrated us- 
ing the cosmic ray signals [23] , as explained in the Ap- 
pendix) , we exclude the presence of GW short signal and 
we establish the upper limit to the amount of isotropic 
GW energy Eqw emitted during the exceptional peri- 
ods. The second algorithm regards the cross-correlation 
function calculated between the EXPLORER and NAU- 
TILUS data streams, as shown in SecjVj Following previ- 
ous methods |24H27| . cumulative analyses are performed 
on the sample of seventy-three measurements obtaining 
several mismatches of the order of 1% probability or less, 
from the background. The statistical significance is eval- 
uated on the basis of the local probability and binomial 
tests of the resulting loudest measurements. In the final 
discussion, Sec |VI[ we provide a summary of the study 
and further comments. 



II. EXPLORER AND NAUTILUS SENSITIVITY 

We consider the EXPLORER and NAUTILUS oper- 
ating during 2004 [3T], using the filtered data set-up as 
the ROG collaboration released f28]. The main char- 
acteristics of the experimental apparatuses, in terms of 
frequency coverage, are shown in table |T] The frequency 
region where the detectors are sensitive covers a range 
of about one hundred hertz, thanks to which the data 
can be sampled at time intervals of 3.2 ms. The data are 
filtered with a filter matched to delta-like signals for the 
detection of short bursts [2^] , typically of the order of a 
few milliseconds. 

Let x{t) be the filtered output of the detector. For well 
behaved noise due to the thermal motion of the oscilla- 
tors and to the electronic noise of the amplifier, the dis- 
tribution of x{t) is normal with zero mean. Its variance 
represents the effective temperature and is indicated 
with Teff. The distribution of x(t) is 
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During the year 2004 T^jf was 2.0 mK for EXPLORER 
and 1.7 mK for NAUTILUS. 



TABLE I: Resonant frequencies of EXPLORER and NAU- 
TILUS during 2004, instrumental bandwidth A/ and relative 
sensitivities. 



detector 


resonances f[Hz] 


Af[Hz] 


Hfnin [ jj^ ] 


EXPLORER 


904.7, 921.3 


8.7 


3.4- 10-22 


NAUTILUS 


926.3, 941.5 


9.6 


3.1 • 10-22 



The resonant detectors are able to measure the Fourier 
component H of the incoming gravitational radiation 
with adimensional amplitude h{t). The apparatuses sen- 
sitivity (SNR=1) to short burst can be expressed as the 
minimum detectable H [30l |3T] 
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where ks is the Boltzmann constant and /o, L and M 
are the resonance frequency, the length and the mass of 
the bar. In general, if we detect a burst with energy Eg, 
expressed in kelvin units, the corresponding value is 
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for optimal orientation. 

Besides the single bursts, this search is suitable for any 
GW transient which shows a nearly flat Fourier spec- 
trum at the two resonant frequencies of each detector. 
The metric perturbation h(t) can either be a millisecond 
pulse, a signal made by a few millisecond cycles, or a 
damped sinusoid signal. In the hypothesis of a signals 
sweeping in frequency through the detector resonances, 
with small decay times (< 50 ms), the filter maintains 
good detection capability [521 - IM| . For a short signal of 
duration time Tg and frequency fg ~ /o, the spectral 
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amplitude, h — \/ J |iJ(/)p df, can be put 



h^HJn/Tg [ 



III. DATA SELECTION 

During the year 2004 several detectors on spacecrafts 
observed a special activity from SGR 1806-20 [? ] culmi- 
nating on December 27*'*. Before the GF, other consider- 
able events occurred, such as the strong outburst October 
5*'', when more than one hundred short bursts were emit- 
ted in a few minutes, involving more than 3 10^^ erg in 
terms of EM energy release |22j. In the last reference, 
there is an accurate description of the event and relative 
light curves of the initial brighter part. The instruments 
(in the 15-200 keV energy range) indicate two very in- 
tense peak clusters, the first one in the vicinity of the 
trigger at 13:56:49 UT, the second one after about three 
minutes. 
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A. On-source region 

Referring to the cited EM triggers, we consider physi- 
cally interesting the following time periods of the filtered 
data of both EXPLORER and NAUTILUS: 

- a) on October 5*'', 13:56:44 - 13:57:10, 26s (starting 
5s before the EM trigger); 

- b) on October 5*'', 13:59:39 - 14:00:25, 46s during the 
second cluster of bright bursts; 

- c) on December 27*'', 21:30:26.14 - 21:30:27.14, Is 
centered to the GF trigger time [SI [7] • 

The choice of the two long first periods is due to the 
special profile of the light curves, but it is also opportune 
to avoid any uncertainty about the time delay between 
the two emissions that make the time reference definition 
difficult in the correlation study. Nevertheless, relatively 
to the GF, it is possible to restrict the measurement pe- 
riod because the impinging time of a light speed signal 
onto the GW detectors is known [14]. Since the three 
different lengths (26s, 46s, Is) don't allow to employ an 
easy algorithm for a uniform analysis from the statistical 
and physical point of view, we divide the two long peri- 
ods into a number of contiguous intervals, 1-second each 
one, obtaining a comprehensive sample of seventy-three 
independent segments for performing measurements. We 
call the GW detector data in this region on-source data. 
All routine checks are done on the selected period, prin- 
cipally by vetoing with a threshold on T^f / , to insure 
good sensitivity and stationary data. In addition, we 
also control the presence of local noise anomalies, and the 
effectiveness of the division into 1-second segments, eval- 
uating the autocorrelation function to the EXPLORER 
and NAUTILUS data of the on-source a and b periods, 
as explained in the Appendix. 



B. The background data 
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FIG. 1: Distribution of the 73 AS parameter values (the 
darkest area), as defined in Sec IV The lightest part is the 



background distribution of the same parameter evaluated on 
the off-source regions. The Kolmogoroff compatibility test 
gives 45% probability. The arrows indicates the AS evaluation 
on the data streams at the time of the GF. 



on source segments, 
indicated with 



The averaged signal parameter is 



AS : 



2n 



K] 



(5) 



For estimating the expected result under the null- 
hypothesis, we use real data but considering time regions 
distant a few hours away the EM triggers. By assigning 
3000 random time steps, we obtain as many off-source Is 
segments, for the first GW detector. To avoid simultane- 
ous signals of any nature, we apply a time shift of several 
hours on the data sequence of second GW detector with 
the respect the first one, thereby obtaining as many as 
3000 off-source segments. All periods are selected vetoing 
by a threshold on Tg^ff, consistently with the on-source 
region selection. 

IV. THE AVERAGE ALGORITHM 

The first algorithm for scanning the on-source data is 
designed to detect an enhancement due to a GW burst. 
Due to the oscillating character of the response, we aver- 
age over Is the amplitude modulus of each sample, Xi for 
EXPLORER and y, for NAUTILUS, belonging to the 



where n — 313, the number of 3.2 ms samples in each 
second. Applying the eq.(|5]) to the 73 pairs on-source 
segments, we obtain the distribution in fig[l] The corre- 
sponding background, the lightest area of fig[l] is cal- 
culated on the 3000x2 segments that have previously 
defined. Cumulative calculations give: ASj^ — 3.29 • 
10^2 vgj.y (,iogg to the AStck value (3.28 • 10"^ 

with root-mean-squared, RMS — 2.8 • 10""^ VK)- Com- 
paring the distributions by the Kolmogoroff test, we ob- 
tain probability 45% that the AS measured distribution 
be similar to the background one. The high compatibil- 
ity value and, as we can easily note, the absence of loud 
events in the distribution, clearly indicate that no ampli- 
tude excess is present on the physical region. Thus, we 
don't retain to perform further statistical investigation 
on this result. Rather, considering the particular phe- 
nomenon and the experimental calibration by the CR ex- 
citation, this measure assumes an important significance 
in terms of upper limit for energy evaluation of emitted 
gravitational radiation. 
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A. Upper limit 

In the present study, two different levels are present in 
terms of source type. 1) Several peaks belonging to the 
outburst on October 5*^, when the EM energy release of 
3 ■ lO-^^ ej.g ^jjg involved. 2) The GF with of 1.6 • lO^^ 
erg of isotropic Eem in the main peak [7j. There are 
four order of magnitude in terms of emitted energy be- 
tween them, so we need to distinguish two measurement 
levels. In Tab. II, corresponding to each category, we 
set out the AS measured parameter. By the cosmic ray 
calibration procedure, as explained in appendix, the in- 
put signal energy values are extracted. From the eqsjS] 
and [4j assuming an incoming burst with duration time 
Tg and / ~ /o, the spectral amplitudes h are evaluated. 
Under the frequentist approach the corresponding 
upper limits are calculated, at 90% confidence level, set- 
ting h = 2.9 ■ 10^21 \/l/Hz, for the period 1), and 

= 4.3 • 10~2i y/l/Hz, for the period 2). 

Using the classical resonant detector cross-section, im- 
der the hypothesis of isotropic emission, the total energy 
carried by the gravitational waves is given by 



E, 



GW 



8G ML^pT, 



37.4 • 10 



49 



1 ms 



' 10 kpc' 



'■Es(mK) [erg] (6) 



where r is the source distance, Tg is the signal duration 
and Eg expressed in mK units. For the categories in 
Tabln) under the hypotheses of an isotropic emission. 



the corresponding Eqw upper limits are 1.8 • 10**^ erg 
and 4 • 10*^ erg, with duration time Tg = 30ms and 
r = 10 kpc. In order to compare the results to previ- 
ous similar measurements, we need to consider the ratio 
•^^"^ The importance of this parameter is well 
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pointed in a series of LIGO analyses regarding the asso- 
ciation with SGRs [TBHin]- Their results are distributed 
over a wide interval comprising almost eight orders of 
magnitude in terms of 7, depending on several signal 
waveform hypotheses and on relative sensitivities. In our 
measurements, 7^"'^° = 6 • 10^, for October outburst, and 
^90% ^ 2.5 • 10^, for the December GF. These values 
are included in the ranges established by previous cited 
measurements . 
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FIG. 2: Distribution of Vmax, the maximum value of the 
cross-correlation function for r = ±16 ms, evaluated on the 
seventy-three on-source regions (the darkest area). The light- 
est part is the distribution of the same parameter evaluated 
on the background 3000 ofT-source regions. The Kolmogoroff 
compatibility test gives 1.3% probability. The arrows indi- 
cates the value at the time of GF. 



TABLE 11: Upper limits for gravitational waves with the EX- 
PLORER and NAUTILUS detectors, using the AS quantity. 



period 
2004 day lenght(s) 


AS 
[VK] 


1^90% 


p90% 
J^GW 

[erg] 


EM 
[erg] 


1) Oct S**" 26-h46 


3.29 ■ 10"^ 


2.9- 10"^^ 


1.8- 10*^ 


3 • 10*2 




±0.3-10"^ 








2) Dec 27*'' 1 


3.11 • 10"^ 


4.3- 10^^^ 


4 • 10*'' 


1.6- 10*'^ 




±0.28.10~^ 









where the filtered data of EXPLORER are indicated with 
x{t), and the filtered data of NAUTILUS, with y{t). x 
is the average value of x, y that of y and e[...] is the 
expectation. We consider the maximum value of the cor- 
relation function 



V. THE CROSS-CORRELATION RESPONSE 

We correlate the data of EXPLORER and NAUTILUS 
processed with the filter matched to delta-like signals by 
taking simultaneous time periods of data of both detec- 
tors lasting one second each (on-source regions) calculat- 
ing the cross-correlation function 

^ £[{x{t + t) - x){y{t ^ t) - y)] 
^e[{x{t)-xYeMt)-m 



rjnax = rnax r(T) r e [—16ms, 16ms] (8) 

In this way we take care of possible small time mis- 
matches between the two detectors and also of the fact 
that the real excitation may not be properly shaped with 
the delta-function used in the matched filter. (In the Ap- 
pendix, we show how this procedure is applied to data 
which include cosmic ray excitations.) The distribution 
of the seventy-three r^ax values is shown in figj2]with the 
relative background. Also for rmax the expected distri- 
bution under the null-hypothesis is empirically evaluated 
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FIG. 3: The cumulative local probability p distribution for 
rmax parameter. The seventy-three on-source measurements 
are represented by the bullets. The continuos line is the null- 
hypothesis Vmax distribution normalized to the number of 
measurements. The dashed line gives the excess needed for 
a 1% confidence level in the null hypothesis. The up-arrow 
indicates the local probability level of rmax at the time of the 
GF. The down-arrow points to the best real result in terms 
of binomial probability P>i8(0.075) = 5.5 ■ lO"*^ 

using the 3000 pair stretches built as described in Sec. 
Ill B. The comparison by Kolmogoroff test gives 1.3% of 
probability that the two distributions be similar. 

A. Statistical tests 

Considering that GW signal from individual Is time 
periods could be likely weak, we test for a cumulative sig- 
nature associated with the whole segment sample. In the 
past, this method has been applied to study correlations 
between EXPLORER and NAUTILUS data and gamma 
ray bursts (GRBs) ^24l[25]. It has been also used and 
developed in the second, third and fourth LIGO science 
run data analysis |26| , and more recently, in the enlarged 
LIGO- Virgo collaboration for GRB search [57]. In the 
last cited papers, we find especially useful the probability 
distribution study and the binomial test, in order to bet- 
ter appraise the statistical significance of the loudest cor- 
relations and to select candidate events. In our case, the 
same method can be helpful to mark time segments that 
particularly deviate from the background. The cumula- 
tive local probability distribution (figj3]) is associated to 
the rmax measurements. The tail of the p local distribu- 
tion corresponds to the loudest rmax values. In the same 
figure the expected distribution under the null hypothesis 
is indicated by the solid line. As easily noted, the p tail 



shows an extended deviation from the background. Cor- 
responding to the measurements with local probability 
p < 0.1, we calculate the binomial probabilities P>i{pi) 
of getting i or more events with p < Pi, obtaining a spec- 
trum which ranges down to 5.5 • 10~^ (at pis — 0.075). 
The measurement relative to the GF segment gives the 
second largest rmax value, setting P2 = 7 ■ 10"'^, and 
P>2 = 0.093. We perform further statistical investiga- 
tions taking into account unforeseeable non-stationary 
effects in the data, and the empirical form of the em- 
ployed algorithm. Sorting k = 20000 different arrange- 
ments, composing each one with 73 random draws from 
the null-hypothesis distribution of figj2] we estimate the 
binomial value P>i{pi,k), and the minimum -P(fc)„i„, 
between the seventy-three values of the fc*'* set. Com- 
paring it with the real best result, we count how many 
times P{k)min ^ 5.5 • 10^^. It happens two times out of 
the 20000 different arrangements. We also calculate the 
expected number of events for having an 1% confidence 
excess with respect the null hypothesis, at given prob- 
abilities p < 0.1. These values are shown in figj3] with 
the dashed line. From the comparison with the rmax on- 
source distribution (the bullets in the same figure), we 
note a discrepancy with the null-hypothesis with proba- 
bilities well below 1%, in the region 0.02 < p < 0.1. 



VI. DISCUSSION 

We studied the EXPLORER and NAUTILUS data in 
correlation with the astrophysical gamma-ray bursts dur- 
ing the year 2004, when an exceptional activity was ob- 
served in terms of electromagnetic emission. In partic- 
ular, we examined the long outburst that occurred on 
October, considering the two time intervals containing 
the major peaks, (respectively of twenty-six and forty- 
six seconds), and Is time interval centered on the out- 
standing flare on December 27*''. On the resulting sam- 
ple of seventy-three segments in total. Is each one, we 
employed two algorithms. The first one consists in aver- 
aging the absolute values of the EXPLORER and NAU- 
TILUS responses on the on source region, and calibrating 
with the transfer function of real delta-signals due to the 
cosmic ray excitations. Assuming the hypothesis of GW 
short pulses, no significant excess is found either with the 
overall statistic sample, nor analyzing the segments with 
loudest amplitude values. We obtain the upper limit ex- 
pressed as amplitude signal or GW energy released by the 
astrophysical source with respect to two energy ranges. 
Considering the bright outburst of October, involving a 
release of ~ 10^^ g^g^ ^.^^^i^ ^90% ^ ^ § . ;^q49 ^^^^ 

With regard to the GF, when an isotropic electromag- 
netic energy of 1.6 10^^ erg is assumed, the upper limit 
is 4 • 10^^ erg. They are comparable with previous values 
obtained by LIGO, in the frequency band 100-1000 Hz. 

We also processed the EXPLORER and NAUTILUS 
data, on the same on-source segments, using the cross- 
correlation algorithm and applying a series of statistical 
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tests. This analysis shows a compatibihty of the order 
of 1% with respect to the nuU hypothesis, for many on 
source segments. Considering the 7 factor, the param- 
eter probing into the GW emission efficiency, we obtain 
~ 10^ and ^ 10'^, implying that the measure results are 
yet far away from the expected level under theoretical 
hypotheses. Therefore, the use of statistical arguments 
alone is not sufficient for deducing a systematic associ- 
ation between the two phenomena, but the study can 
turn out useful for further investigations. More obser- 
vations are necessary in the direction of these intense 
sources. From this point of view, the recent improve- 
ments in gamma-astronomy are very encouraging. Over 
last few years, various other important detectors have 
been launched, such as the AGILE mission [3B and Fermi 
Gamma-ray Space Telescope |37j . Other magnetars were 
also detected, such as the SGR 0501-^4516 [38], the SGR 
0418-^5729 [391, and the very recent SGR 1833-0832 [40]. 
Thus, a fair amount of solid evidence exists for a more 
definite scientific strategy, studying multi-messenger phe- 
nomena. 
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FIG. 4: Autocorrelation functions for both EXPLORER (up) 
and NAUTILUS (down) data. The checks are performed us- 
ing fifty seconds overlapping the on-source periods a and 6 
(respectively on the left and on the right of the figure) both 
belonging to the October 5*'' outburst. 



VIII. APPENDIX 



Autocorrelation test 



Performing the GW data analysis, three time inter- 
vals are chosen, as explained in the Sec |III| The first two 
periods are a few ten seconds long each one, but they 
have been divided into 1-second segments for the anal- 
ysis. To check the presence of local noise structures, as 
little glitch sequences or long bumps, that could invali- 
date the measurement also from the statistical point of 
view, we evaluated the autocorrelation function to the 
EXPLORER and NAUTILUS data of the on-source a 
and b periods chosen in the Sec |III[ The result is shown 
in the fig |VIII C Although the checks are performed us- 
ing fifty seconds of data for each period, the graphics are 
presented up to one second of time shift r, given that 
the autocorrelation function drops within ~ 200 ms , as 
we easily note, and the same behaviors continue on the 
whole periods. 



B. Cosmic rays detection 

The test validity is performed employing the real sig- 
nals generated by extensive air showers (EASs) [23], and 
the fihered data of EXPLORER and NAUTILUS, all reL 
ative to the year 2004 [21, ,28^ . Two events were selected. 



each for the single detector, with highest interaction en- 
ergy with relative bar, in order to be independent of the 
noise. The responses of the two detectors are shown in 
figj5j The selected EASs have respectively 3561 and 3128 
of particle multiplicity measured in the lower part of the 
cosmic ray detectors. The vibrational energy peaks are 
correspondently 2.47Ar and 1.08Ar, several hundred times 
the background, indicating a total energy of the order of 
10 TeV deposited in each bar [23] . 



C. Average algorithm calibration 

We use the typical envelope (see figjs]) for simulating 
delta-like signals. Obviously, we opportunely reduce each 
signal by scale factor depending on the choice of the am- 
plitude signal. Then, we add the constructed signal to 
the centre time of 300x2 off-source regions using 10% 
of the background employed in the Sec jlll B] In terms 
of energy amplitude six values are chosen in the range 
0.5 — 30 mK, and for each level, the average algorithm 
is applied obtaining the corresponding AS responses ap- 
plying the eq[5j The calibration curve is reported in the 
figj6j At Eg — 0, the ASbck value is reported, as calcu- 
lated in SecHvl 
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D. Cross-correlation response 
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FIG. 5: Two selected cosmic ray showers producing in the 
two detectors signals much larger than noise. On the ordinate 
axis we show the ^/~E, with the energy E in kelvin units. We 
notice the different resonance frequencies in the two detectors. 
The time of the two signals are shifted to an arbitrary zero, 
in order to simulate two simultaneous signals. 
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FIG. 6: The AS response, in squared kelvin units, versus the 
energy of cosmic ray calibration signals. Fore each point, 300 
off-source measurements are been employed and averaged. At 
Es = 0, the ASbck value is reported. The continuos line is a 
polynomial fit. 



To study the cross-correlation response, we use the 
EASs selected before and the corresponding GW data 
stretches showed in fig[5l The relative cross-correlation 
function is shown in fig[7| We note the oscillating be- 
haviour of the maximum response envelope in the tem- 
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FIG. 7: The correlation function of the signals shown in the 
previous fig[5] versus the time shift. The zero on the abscissa 
is arbitrary. 



poral shift window of a few ten-milliseconds at zero time 
shift. This can compromise the efhcicncy of the cross- 
correlaction detection algorithm, especially in very low 
signal-to noise-ratio conditions. To counter this, we 
adopt the rmax parameter, the maximum modulus of 
the cross-correlation function in the interval of ±16 ms 
around the zero shift. 
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